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We present the results of our calculation which has been performed to study the nuclear effects 
' in the quasielastic, inelastic and deep inelastic scattering of neutrinos(antineutrinos) from nuclear 

targets. These calculations are done in the local density approximation. We take into account the 
effect of Pauli blocking, Fermi motion, Coulomb effect, renormalization of weak transition strengths 
in the nuclear medium in the case of the quasielastic reaction. The inelastic reaction leading to 
production of pions is calculated in a A- dominance model taking into account the renormalization 
of A properties in the nuclear medium and the final state interaction effects of the outgoing pions 
with the residual nucleus. We discuss the nuclear effects in the F$~(x) structure function in the 
deep inelastic neutrino(antineutrino) reaction using a relativistic framework to describe the nucleon 
q . spectral function in the nucleus. 

£H 1 PACS numbers: 12.15.-y,13.15.+g,13.60.Rj,23.40.Bw,25.30.Pt 

It is now well established that neutrinos oscillate. The next target of the experimentalist is to determine precisely the 
various parameters of neutrino mass matrix given by Pontecarvo-Maki-Nakagawa-Sakata (PMNS), absolute masses 
of different flavors of neutrinos, pattern of the known neutrino mass differences, CP violation in neutrino sector, 
^2 ' e tc- To address some of these problems several experiments like CNGS, MINOS and SciBooNE are taking data. 
The eperimental analyses of neutrino oscillation data are going on, for example, at K2K and MiniBooNE, and 
several experiments are planned to be done in future like T2K and NOi^A. Besides the accelerator experiments, 
experiments with neutrinos from v- factories, /3-beams, etc., are also planned, as well as some experiments with 
natural j/-sources like solar neutrinos, atmospheric neutrinos, or (anti)neutrinos from nuclear reactors are also planned. 
These experiments use various nuclear targets like 12 C by MiniBooNE, SciBooNE, MINERVA and NO^A, le O by 
SuperKamiokande, T2K, UNO, Hyper-K, K2K and MEMPHYS, 40 Ar by ICARUS and NOz^A, 56 Fe by MINOS, INO, 
MINERVA, and 208 Pb by CNGS, MINERVA and OPERA collaborations. Most of these experiments are being done in 
the neutrino energy region of E v < 2GeV, for example at MiniBooNE, the average energy (< E v >) is 750 MeV while 
at K2K it is 1.3 GeV. At the energies of a few GeV the contribution to the cross section comes from the quasielastic, 
inelastic and deep inelastic processes and the analysis of the data of neutrino experiments is based on Monte Carlo 
generator of events like NUANCE, NEUGEN, NEUT, etc. In these Monte Carlo generators, neutrino cross sections 
are used which are based on the model of Llewellyn Smith pland Smith and Moniz 0] for the quasielastic reactions, 
Rein and Sehgal [3;] for the inelastic reactions and GRV98 4] along with the modifications suggested by Bodek et 
al. [j| for the deep inelastic reactions. 

The importance of a better knowledge of neutrino-nucleus cross section to be used in the Monte Carlo generators 
has been realized and discussed in a series of neutrino conferences like Nulnt, NuFact, NOW, etc. There are, now, 
various theoretical calculations for the quasielastic process which make use of nuclear models like shell model with 
pairing correlations, random phase approximation, relativistic mean field approximation, etc. In the case of inelastic 
neutrino-nuclear reactions pion production processes have been studied. These are generally studied in a A-dominance 
model in which pions are dominantly produced through the excitation of A and its subsequent decay leading to pions. 
Some of the calculations have also been done by taking background as well as higher resonance terms. In the case 
of deep inelastic scattering of neutrinos (antincutrinos) from nuclear targets, there are very few calculations where 
the dynamical origin of the nuclear medium effects have been studied. In some phenomcnological analyses, nuclear 
medium effects have been described in terms of few parameters which arc determined from fitting the experimental 
data of charged leptons and neutrino (antineutrino) deep inelastic scattering from various nuclear targets. 

Here we discuss the nuclear effects in the quasielastic reaction, inelastic one 7r production in the A dominance 
model and the nuclear effects on the F-^(x,Q 2 ) structure function in the deep inelastic reaction. These calculations 
are done in the local density approximation. For the quasielastic process this model takes into account the effect of 
Pauli blocking, Fermi motion, Coulomb effect, renormalization of weak transition strengths in the nuclear medium. 
The inelastic reaction leading to production of pions is calculated in a A- dominance model taking into account the 
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renormalization of A properties in the nuclear medium and the final state interaction effects on the outgoing pions. 
We discuss the nuclear effects in the F^(x) structure function in the deep inelastic neutrino(antineutrino) reaction 
using a relativistic framework to describe the nucleon spectral function in the nucleus. The details of these calculations 
may be found in the Refs. @, 0] for the quasielastic process, Refs. [j| Q for the inelastic process and Ref. [l(| for 
the deep inelastic process. Similar calculations for nuclear effects in the quasielastic and inelastic processes have also 
been recently done by many other groups [lll|-(l7j. In the following we describe, in brief, the formalism for calculating 
the nuclear effects in quasielastic, inelastic and deep inelastic processes in Sec.I, and present the numerical results in 
Sec. 2 with concluding remarks given in Sec. 3. 



I. NEUTRINO NUCLEUS REACTIONS 



A. QUASIELASTIC REACTION 

The basic reaction for the quasielastic process is a neutrino interacting with a neutron inside the nucleus which is 
given by 

!/„(*;)+ n(p)-» ti-(k')+p(p') (1) 

The cross section for quasi-elastic charged lepton production is calculated in the local density approximation by taking 
into account the Fermi motion and the Pauli blocking effects through the imaginary part of the Lindhard function 
for the particle hole excitations in the nuclear medium. The renormalization of the weak transition strengths are 
calculated in the random phase approximation(RPA) through the interaction of the p-h excitations as they propagate 
in the nuclear medium using a nucleon-nucleon potential described by pion and rho exchanges. The effect of the 
Coulomb distortion of muon in the field of final nucleus is taken into account using a local version of the modified 
effective momentum approximation. 

The total cross section o(E v ) for the charged current neutrino induced reaction on a nucleon inside the nucleus in 
a local Fermi gas model is written as p]]: 

a ( Ev ) = - 2Gf2c ° s2 ° c f ma % 2 dr f P " p, 2 dp, C d(cos9)-^—L tlu J^ImU N (q 0) q). (2) 

7T Jr min Jp™'» J -I ^u^il 

where = £ L„Lj and J»» = £ £ J»J"*. 

The leptonic current L M and the hadronic current J M are given by 

L„ = u(k'h„{l - j 5 )u{k) (3) 



= faV + F 2 (q 2 )ia^^ + F A (q 2 )n 5 + F P (q 2 )q^ 5 ]u(p). (4) 

where q(= k — k') is the four momentum transfer, M is the mass of the nucleon, Gf{= 1.16637 x 10 _5 GeT^ -2 ) is the 
Fermi coupling constant and 9 is the lepton angle. Un is the Lindhard function for the particle hole excitation 6}. The 
form factors F% , F 2l Fa and Fp are isovector electroweak form factors and for our numerical calculations we have used 
the parameterisation of Bradford et al. [l|| with axial dipole mass Myt=1.05GeV and vector dipole mass My =0.84GeV. 
Inside the nucleus, the Q-value of the reaction and Coulomb distortion of outgoing lepton are taken into account by 
modifying the imaginary part of the Lindhard function 7m[/jv(9oi Q) by ImU^{q^ — V c (r) — Q,q). Furthermore, 
the renormalization of weak transition strength in the nuclear medium in a random phase approximation(RPA) is 
taken into account by considering the propagation of particle hole(ph) as well as delta-hole(A/i) excitations. These 
considerations lead to modified terms involving the bilinear terms in the weak coupling constant in the hadronic tensor 
Jrpa f° r wn ich expressions are given in Ref. [2|. 



B. INELASTIC RESONANCE PRODUCTION OF PIONS 



The basic reaction for the inelastic one pion production in nuclei, for a neutrino interacting with a nucleon inside 
a nuclear target is given by 

v^k) + N(p) -> p-{k r ) + N'{p') + it + {K) N,N'=p/n (5) 
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The cross sections for pion production is calculated using the A dominance model. In this model, the weak hadronic 
currents interacting with the nucleons in the nuclear medium excite a A resonance which decays into pions and 
nucleons. The pions interact with the nucleus inside the nuclear medium before coming out. The final state interaction 
of pions leading to elastic, charge exchange scattering and the absorption of pions lead to reduction of pion yield. The 
nuclear medium effects on A properties lead to modification in its mass and width which have been discussed earlier 
by Oset et al. to explain the pion and electron induced pion production processes from nuclei. 

In the local density approximation the expression for the total cross section for the charged current one pion 
production is given by 



i , _ 1 



( 47r ) 5 Jr min Kr " VJ ' 9' 



a = 7Z^5 / (Pp( r ) + ^Pn(r))df I ^ dQ 2 J dk' J d{cosO n ) J d(f> n X 

-EEi-^i 2 (6) 



MEZE t E , p + En (i_WL cos{ 

where the proton density p p (r) = ^p(r) and the neutron density p n {f) = ~2 Z p(t) with p(r) as the nuclear density 
taken as 3-parameter Fermi Density taken from Ref . [2l| . The transition matrix element M. ft is given by 

M fi = V3 Gf * LNA *(P)KV aX Xa L a u(p) (7) 

V 2 m ir 

where L a is the leptonic current defined by Eq.Q, a = cos9 c and 0? a = O v a + P A a for the charged current 
induced tt^ production process while for the neutral current induced tt° production process a=l and 0@ a = (1 — 
2sin 2 9w)Oy a + C^ Q - O v a and are the vector and axial vector N-A transition operators given by 

°y = i^w 1 ^ A - Aq) + ^ 1 ^ aPq ■ p ~ qPpa) + ■ p ~ ^ ^ 

and 

< = ~ «V) + <tf («V + (9) 

where (q 2 ) and Cf(q 2 ) are the vector and axial vector transition form factors and for our numerical calculations 
these have been taken from the work of Lalakulich et al. [20| . Ow is the weak mixing angle. V iyX is the A propagator 
in momentum space given by 

pcrA 

V aX = (10) 

P 2 - M 2 A + iM A T 1 ' 

where P is the spin-3/2 projection operator given by 

_ E ^ . ( /> + A/ A ) (p- - + l^jjp^ - \f+) (ID 

spins \ £i s 

and the delta decay width T is taken to be an energy dependent P-wave decay width taken as [l9| : 

1 ( IttNa\ 2 M,_ 3 



T(W) = - ^ -l|q cm | J (12) 
bir \ J W 

1 9cm | is the pion momentum in the rest frame of the resonance and W is the center of mass energy. 

Inside the nuclear medium the mass and width of delta are modified which in the present calculation are taken 
into account by using a modified mass — > + ReY>& and modified width r a — > ~ 2Pm£A from the model 
developed by Oset et al. [li| . where f a is reduced width of A due to Pauli blocking of nucleons in the A — > Ntt 
decay and Sa is the self energy of A calculated in nuclear many body theory using local density approximation. The 
expressions of Pe^A and JtoSa are taken from the Ref. (l9| . The pions produced in this process are scattered and 
absorbed in the nuclear medium. This is treated in a Monte Carlo simulation which has been taken from the Ref. 12211. 
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C. DEEP INELASTIC REACTION 

The basic process for a neutrino interacting with a nucleon inside the nucleus is 

u ll {k)+N{p)^ l x-{k')+X{p'). (13) 

where X is the jet of partons. 

The differential scattering cross section for the deep inelastic scattering of (anti)ncutrinos from unpolarized nucleons 
in the limit of lepton mass mi — > 0, is described in terms of three structure functions, Fi(x,Q 2 ), F^ix.Q 2 ) and 
F 3 (x,Q 2 ), where x = ^7 = — ^ is the Bjorken variable, v and q being the energy and momentum transfer of 
leptons. In the asymptotic region of Bjorken scaling i.e. Q 2 — > 00, v — > 00, x finite, all the structure functions depend 
only on the Bjorken variable x. In this scaling limit, F"(x) and F%(x) are related by the Callan-Gross relation [23| 
leading to only two independent structure functions F^{x) and F 3 (x) which are determined from the experimental 
data on deep inelastic scattering of (anti) neutrinos in the asymptotic region. 

We have studied nuclear medium effects on the nucleon structure function F 3 A {x,Q 2 ) in iron using spectral function 
to describe the momentum distribution of nucleons in the nucleus. The spectral function has been calculated using the 
Lehmann's representation for the relativistic nucleon propagator and nuclear many body theory is used to calculate 
it for an interacting Fermi sea in nuclear matter. A local density approximation is then applied to translate these 
results to finite nuclei [24|. Here we consider the modifications of nucleonic contributions to F 3 A (x,Q 2 ) arising due to 
binding energy, off mass shell and Fermi motion of the nucleon in the nuclear medium which dominate in the region of 
x >0.3. In this model, in the region of 0.3> x >0.1, corresponding to the anti-shadowing region, the nuclear medium 
modification effects on F 3 A (x,Q 2 ) are expected to be small due to vanishing of the pion contribution and we have 
not considered the shadowing region of 0.0< x <0.1. Therefore, our results should be able to describe the dominant 
contribution of nuclear medium effects to F 3 {x,Q 2 ) in the range of 0.1< x <1. 

The average structure function F 3 (x) on isoscalar nucleon target defined as 

F 3 N (x)= 1 ~(F 3 »»+Fr) 

is given by F$ (x) = \u v (x) + d v (x) + s(x) — s(x) + c(x) — c(x)}, where u v {x) and d v (x) are the valence quark parton 
distributions. For an isoscalar target and a symmetric sea, F 3 N (x) structure function is given in terms of valence 
quarks u v and d v which satisfy the Gross-Llewellyn Smith sum rule (25j: 

f F 3 N {x)dx = 3. (14) 
Jo 

In the local density approximation the reaction given by Eq. (|13p takes place at a point r, lying inside the nucleus in 
a volume element d 3 r with local density p p (r) and p n (r) corresponding to the proton and neutron. The expression 
for F$(x, Q 2 ) in the nuclear medium is given by [lfj| |: 

F 3 A (x,Q 2 ) = Aj d?r J -0, J 11 dp S h (p°,p,p(r))F(p,Q 2 )F 3 N (x N ,Q 2 ), (15) 
where xn is the Bjorken variable expressed in terms of the nucleon variables, (p° , p), in the nucleus i.e. xn = 

= £(p) { (^^h ) 5 7 = I 1 + ~W J and 

Sh(u, p) is the hole spectral function, the expression for which is taken from Ref. (26j . 

II. RESULTS AND DISCUSSION 

A. QUASIELASTIC REACTION 

In Fig.l, we present the ratio R of the charged current quasielastic lepton production cross section to the cross 
section on free nucleon defined as R = CT (y> ee ) as a function of neutrino(Fig.la) and antineutrino(Fig.lb) energies 
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FIG. 1: R = ffeS vs Neutrino (Antineutrino) Energy. 
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FIG. 2: -j^- vs p M for the f M (f M ) induced reactions on 12 C target at E v — lGeV. 



when cr( 12 C) is calculated using Eq.(j2]). We find that with the incorporation of nuclear medium effects without the 
RPA correlations the reduction in the cross section is around 45% at i?„=0.2GeV, 16% at E v =0AGeV, 10% at E v w 
l.OGeV and 7% at i?„=2GeV from the cross sections calculated for the free case. However, when we also incorporate 
the RPA effects, the total reduction in the cross section is around 40% at i?„=0.4GeV, 20% at 2£„=1.0GeV, 18% at 
E„=2GeV from the cross sections calculated for the free case. In the case of antineutrinos these reductions are larger 
as shown in Fig. 1 (b) . We have compared our results with the results obtained in the Fermi gas model which has been 
used in the NUANCE Monte Carlo generator [27[ by the MiniBooNE collaboration [11]. We find that the present 
results in the local Fermi gas model are similar to the results used in the NUANCE generator, but when RPA effects 
are included the cross sections are reduced. 

In Fig. 2, we have shown the results for the lepton momentum distribution for the and induced charged 
current quasielastic processes. We find that when nuclear medium effects are taken into account there is a reduction 
as well as shift in the peak region towards the lower value of lepton momentum. This reduction in when calculated 
in the local Fermi gas model without the RPA correlation effects as compared to the cross section calculated without 
the nuclear medium effects is around 10% in the peak region of lepton momentum, which further reduces by around 
30% when RPA effects are also taken into account. In the case of antineutrino the reduction in 4 s - in the local Fermi 

gas model is around 30% which further reduces by 30% when RPA effects are also taken into account. 
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FIG. 3: -j^j and for the z-vt^p) induced charged current one n (tt ) process on C target at E v = lGeV 



B. INELASTIC REACTION 

In Fig. 3, we present the results for Q 2 -distribution and momentum distribution for the charged current 

v tiiyti) induced one 7r + (7r~) production cross section. These results are presented for the differential scattering cross 
section calculated with and without the nuclear medium effects and with nuclear medium effects including the pion 
absorption effects. For the Q 2 - distribution shown in Fig. 3a, we find that the reduction in the cross section as 
compared to the cross section calculated without the nuclear medium effects is around 35% in the peak region. When 
pion absorpion effects are also taken into account there is a further reduction of around 15%. The results for the 
antineutrino induced one 7r~ production cross section are quantitatively similar in nature but quantitatively we find 
that the peak shifts towards a slightly lower value of Q 2 . In Fig. 3b, the results for the pion momentum distribution 
have been shown. We find that in the peak region the reduction in the cross section is around 40% when nuclear 
medium effects arc taken into account, which further reduces by about 15% when pion absorption effects are also taken 
into account. In Fig. 4, we present the results for the neutral current v{v) induced one tt° production cross section. 
These results are presented for the pion momentum and angular distributions dc ^ g , with and without the nuclear 
medium effects and with nuclear medium and pion absorption effects. For the pion momentum distribution shown in 
Fig. 4a, we find that in the peak region the reduction in the cross section is around 40% when nuclear medium effects 
are taken into account, which further reduces by about 15% when pion absorption effects are also taken into account. 
The results with antineutrinos are similar in nature, except that in the case of D, the angular distribution are more 



7 



£ 

o 

t 2 
o 



Oh 

-a 

C 1 







~i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — r 

Neutrino 



Free 

Medium Effects(ME) 
ME + Pion Abspn 




l — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — r 



Antineutrino 



Free 

Medium Effects(ME) 
ME + Pion Abspn 




200 400 600 800 1000 u 200 400 600 800 1000 

p^MeV) p (MeV) 



1.6 



B 1.2 

o 



O 

o 

"B 0.4 

-a 







_ 1 1 


1 1 1 1 1 1 1 1 1 


i I i i i , i 

* ♦ _ 




Neutrino 






■ ■ ■ Free 


• 




— ■ Medium Effects(ME) 


* _ 




— ME + Pion Abspn 


; ' \ 




(b) 














♦ *s 










• 


Syr 








1 1 


i i i i i i i i i 


i l i i i l i 



0.6 



-0.8 -0.4 0.4 O.i 

cos 6 



0.2 







1 


1 1 1 1 1 1 1 1 1 1 


1 1 | 1 1 1 | 1 




Antineutrino 






■ ■ ■ Free 


» 




— ■ Medium Effects(ME) 






— ME + Pion Abspn 


•* n 




(b) 








/ '/ 






'/ - 
















...•*' 




1 


i i i i i i i i i i 


1 1 1 1 1 1 1 1 




-0.8 -0.4 


0.4 0.8 



cos 6 



FIG. 4: -j 2 - and — for the neutal current neutrino(antineutrino) induced n° production on 12 C target at E v = lGeV. 



forward peaked than in the case of v. 

In Fig. 5, we present the results for the total scattering cross section a for charged current v^Va) induced one 
n + (ir~) production cross section. These results have been presented for the cross sections calculated without(with) 
the nuclear medium effects and also when pion absorption effect is included along with the nuclear medium effects. 
We find that with the inclusion of nuclear medium effects the reduction in the cross section from the cross section 
calculated without the nuclear medium effects for neutrino energies between 1-2 GeV is 30-35% which further reduces 
by 15% when pion absorption effects are also taken into account. The results with antineutrinos are similar in nature. 



C. DEEP INELASTIC REACTION 

In Fig. 6(a), we compare our results for K(x,Q 2 ) at Q 2 — 5GeV 2 , where R = ^pN^ql^ , with the results of Tzanov 

ct al. 2t| JKulagin and Petti [3(j, Kulagin [3l[ and Hirai et al. [32j]. While the work of Kulagin [3l[ and Kulagin and 
Petti 3fJ,[33| use a nuclear model to calculate the nuclear effects which shows a Q 2 dependence, the work of Tzanov et 
al. [2£ | and Hirai et al. [13] are phenomenological analyses, which assume the nuclear effects to be independent of Q 2 . 
We find a suppression in F^(x,Q 2 ) for x <0.7 and an enhancement thereafter, which are respectively smaller than 
the results of Kulagin [3l| , but are larger than the recent results of Kulagin and Petti [33[ . It should be noted that 
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FIG. 5: a for (u^) induced charged current one ty + (ty ) production on 12 C target. 




FIG. 6: (a)Results for the ratio R = fS.N^\ at Q 2 ~ 5GeV 2 by different authors. Solid line: Result of the present calculation 

using MRST2004 NNLO parton distribution function; short dashed line: NuTeV Collaboration [29( ] . dashed-double dotted line 
is the result of Kulagin [3lj], double dashed-dotted line: Hirai et al. [3^ ] and the results of Kulagin and Petti are shown by 
dashed-dotted line.(b)AGLS = ±(3 - J,, 1 Fi(x, Q 2 )dx) vs Q 2 . 



these latter results [33[ give suppression in the region of 0.4< x <0.8 and enhancement for x >0.8, which are smaller 
than the present results and the results obtained earlier in Ref. [3l|. When compared with the results of Tzanov et 
al. [29( and Hirai et al. [321 ]. we find a smaller suppression in the region 0.5< x <0.7. In the region 0.7< x <0.8, we 
find an enhancement while they obtain a suppression. In Fig. 6(b), we show the Q 2 dependence of the nuclear effects 
of the GLS integral, where we plot AGLS=|(3-/ D ] F£(x,Q 2 )dx) function of Q 2 . The experimental results from 
CCFR collaborations [HI, CHARM collaborations [35[ and IHEP-JINR collaborations [11] are also shown. The Q 2 
behaviour of AGLS has been found to be in reasonable agreement with the present available experimental results. In 
this figure, we have also shown the theoretical results obtained by Qiu and Vitev [37|- Our results are in agreement 



with the results of Qiu and Vitev [37J for Q 2 > 5GeV 2 where theoretically the suppression is found to be larger 



than the experimental results. For Q 2 < 5GeV 2 , we find a larger suppression compared to the central value of the 
experimental result and both theoretical values are within the experimental errors. 
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III. CONCLUSIONS 



We will like to conclude that : 

(i) In the case of charged current quasielastic lepton production, the role of nuclear medium effects like Pauli 
blocking, Fermi motion is to reduce the cross section. When nuclear correlations in the nuclear medium are taken into 
account in a Random Phase Aprroximation (RPA) there is further reduction in the cross section. The total reduction 
in the cross section is about 20% in case of v and slightly larger for v around — lGcV. 

(ii) In the case of charged current one pion production, the nuclear medium and pion absorption effects lead to a 
reduction in the cross section about 45% at E u =\ GeV and give appreciable distortion in the energy and angular 
distribution of pions and leptons. 

(iii) In the case of deep inelastic scattering, the nuclear effects decrease the value of the structure function F%(x, Q 2 ) 
in the iron nucleus for x < x m i n —0.7 and increase it at higher x > x m i n . In general nuclear medium effects decrease 
the value of GLS integral for all Q 2 . 
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